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Abstract 

We study a hard regime of stimulation of two-frequency oscilla- 
tions in the main resonance equation with a fast oscillating external 
force: eitp' + IV'PV' = sxp (it^/(2e)), < e <C 1. This phenomenon 
is caused by resonance between an eigenmode and the external force. 
The asymptotic solution before, inside and after the resonance layer 
is studied in detail and matched. 

1 Introduction 

In this paper we investigate the hard mode of stimulating of two phase os- 
cillations. We study this phenomenon in asymptotic solution of ordinary 
nonlinear differential equation under fast oscillating external force: 

2^ 



/it \ 

eillj' + = exp i^—j 



where < £ ^ 1 - a small parameter. 

The solution of the equation (0) constructed in this paper oscillates with 
the frequency of the perturbation when t > t*. An amplitude of oscillations 
changes slowly. This system changing corresponds to the stimulated oscilla- 
tions of nonlinear equations. The frequency of the oscillations changes with 
time. From the other hand, the eigen-frequency of the nonlinear equation de- 
pends on the amplitude. At certain moment t = the eigen-frequency of the 
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oscillations becomes equal to the external force frequency. A resonance takes 
place in the system. This leads to hard loss of stability of the eigenmode of 
the oscillations and the system gets out of the resonance. The asymptotic 
solution of equation (|l|) becomes two-phase. One phase relies to the oscilla- 
tions stimulated by the perturbation force and another one corresponds to 
the oscillations stimulated by the transition over the resonance. 

To explain the nature of the solution bifurcation and to construct the 
asymptotic solution it is more convenient to investigate the equation for the 
amplitude of the oscillations U = ■j/'exp ( — it^ /{2e)) 

eiU' +\U\'^U -tU = 1. (2) 

If we begin to study the amplitude U and consider the equation (0) then 
constructing of the asymptotic solution is equivalent to investigation of the 
bifurcation of the slowly varied equilibrium of the equation (0). The theory 
of the bifurcation usually investigates the behavior of a solution depending 
on external parameter which is not connected with a variable of a differential 
equation p|, ^. 

The bifurcation of the equilibrium with varying parameter is well known 
in physics 0. This bifurcation is usually illustrated by loss of stability of 
the equilibrium of mechanical system, which is described by the second order 
differential equation with slowly varying coefficients. In this case the hard 
loss of stability takes place. The numeric evaluations give the picture: 




Figure 1: The bifurcation solution of the equation (^. 

A description of the internal asymptotic structure of the solution in the 
case of the hard loss of the stability is connected with a construction of 
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the asymptotics for the solution of the Painleve-1 equation [|]. Besides the 
solution has a complicated asymptotic structure in the transition layer where 
the main term of the asymptotics is defined by four various expansions of 
different types. There are: a special solution of the Painleve-1 equation, a 
sequence of separatrix solutions of the homogeneous equation with "frozen" 
coefficients, a sequence of solutions of the Weierstrass equation with the 
parameter (72 = and a sequence of solutions of the Weierstrass equation 
with the parameter (72 7^ 0. 

Two of the first internal layers were found and studied in [^. Later, in 0, 
[0 the asymptotic structure of the transition layer was studied in detail. It 
was shown that actually it is necessary to study the sequence of the separatrix 
solutions of the homogeneous equation with "frozen" coefficients. The third 
internal layer was found in [0. 

However, these three layers are not enough for a passage through the 
bifurcation interval. There exists a layer which is defined by a solution of 
the Weierstrass equation with the parameter §2 7^ 0. This layer is found 
and studied in this work due to the matching of all asymptotic solutions 
before inside and after the bifurcation layer. The another new result of this 
paper is the calculating of internal variables in a sequence of layers, which are 
connected with the separatrix solutions of the homogeneous equation with 
"frozen" coefficients. For the Painleve-2 equation with a small parameter 
at a derivative these calculations have been made in the work of one of the 
authors of this paper in 

Among other works devoted to the investigation of the asymptotics solu- 
tions for nonlinear differential equations with variable coefficients it is nec- 



essary to note the investigations of passage through a separatrix in p[ -|12 



where there were no confluent of the slowly varying equilibrium. In par- 
ticular, in work |[T2[ the passage through a separatrix of equation (Q) with 
constant coefficients and small dissipation was studied. In the case consid- 
ered in this work the solution of the equation for the amplitude also passes 
through a separatrix, but at more complicated confluent saddle-center equi- 
librium. 

Here, the asymptotic solution for the equation (|l]) in the domain where 
the amplitude varies slowly [t > t*) is constructed by perturbation theory 
method and in domain where the amplitude fast oscillates — by the Krylov- 



Bogolyubov method [0, |14[. Thus, we reproduce the elegant formulas for 
the asymptotic solutions of the equation (|1|) obtained in the work |T^. All 
asymptotics are matched [0. 

The structure of this paper is following. In the second section the prob- 
lem is formulated. In the third one, the obtained results are represented. 
The fourth section is devoted to a construction of an algebraic asymptotic 
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solution for the equation (^. The fifth section includes an investigation of a 
bifurcation layer and studying of the sequences of the internal expansions. In 
the sixth section the fast oscillating solution of the equation (|^) is constructed 
and its asymptotics at the approach to a bifurcation point is written. 

2 Studied problem and its justification 
2.1 Justification 

The general solution of the equation (|l]) generally speaking is unknown. But 
we can construct a set of formal asymptotic solutions on a small parameter 
e. The simplest kind of the solutions is the solutions oscillating with the 
frequency of the external force. To obtain this kind of solutions it is neces- 
sary to proceed the amphtude equation (Q). Further following the natural 
suggestion about boundedness of the derivative in the equation (H) we can 



obtain a nonlinear algebraic equation for the main term of asymptotics U (t): 

\U\^U -tU= 1. (3) 

The number of the roots of this algebraic equation depends on a parameter 
t. The roots can be written explicitly. There exist a value of the parameter 
t equals to = 3(1/2)^/^ so that the equation (|) has three real roots at 
t > t^:. At t = = 3(1/2)^/^ there is one simple root and one multiple root 

= —(1/2)^/'^. At t < the equation (^) has the alone root. 

Let us consider the domain t > where there are three different roots 
of the equation (^). Denote them Uk, and U^it) < U2{t) < Ui{t). These 
roots correspond to slowly varying equilibriums of the equation (Q). Two of 
the equilibriums are stable centers. The formal asymptotic solutions with 
the main terms Ui(t), U2(t) correspond to these stable equilibriums. The 
third equilibrium is a saddle. There exists a parameter value t = t* at which 
one of the centers coalesces with the saddle. At that moment the saddle- 
center bifurcation takes place. At t < there exists only one slowly varying 
equilibrium. 

This bifurcation may by explained on the example of an autonomous 
equation with a "frozen" coefficient T: 

iV + {\V\'^ -T)V = 1. 

On the figure 2 one can see the phase plane of this equation. On the 
left picture one equilibrium exists. On the middle picture one can find two 
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Figure 2: The phase portraits for the equations with the "frozen" coefficient 
at T < t*, T = and T > t*. 

equihbrium positions and at last on the right picture one can find three 
equihbrium positions. 

Choosing Uk as the main term of the asymptotics in the domain t > t^, 
one can construct three different formal asymptotic solutions in the form: 

oo 

U{t-e) = Y.e''U{t). (4) 

n=0 

Each of them corresponds to the one of the slowly varying equilibrium. The 
formal solution with the main term Ui (t) doesn't change the structure when 
the variable t changes and we do not investigate this case. Solution based 
on 1/3(1) is unstable on the whole interval t > and it is also out of our 
consideration. We will investigate the asymptotic solution based on U2(t) and 
the behavior of this asymptotic solution after the saddle-center bifurcation. 

2.2 Statement of the problem 

In this work we solve the problem on constructing of the formal asymptotic 
solution of the equation (|l]) in the interval t G [t* — C, + C] where C = 
const > uniform on e . We suppose that the solution in the domain t > 
has the form 

^(t,e)=exp — e ^ 0, 

„=o 

where U it) = U2{t). 

3 Result 

The combined asymptotic solution of the equation (0) in an interval t G 
[t^ — C,t^ + C] is constructed. The solution has the different asymptotic 
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structure in different parts of tlie interval. Tlie main result of tliis work is 
following: firstly, the domain where constructed combined solution is valid 
completely covers the interval [t*— C, t*+C] and secondly, all asymptotics are 
matched. In this section we only represent the form of asymptotics and leader 
terms. The explicit formulas can be found in the corresponding sections of 
the paper. 

In the domain (t — t^:)e~'^^^ ^ 1 the asymptotic solution has the form: 

^^^^ n=0 

The leader term of the asymptotics is equal to the middle root of algebraic 

n 

equation d): U (t) = U2{t). The correction terms U (t) are algebraic func- 
tions of t. 

In the domain It — -C 1 the asymptotic solution is defined by four 
various expansions of different types. First of them is: 

^P = (u.+e'/' f .^"/^ a ir)+^e'/' £ .^"/^ P (r)) exp (|^), e 0, (5) 

^ n=0 n=0 ^ y^Sj 

where is multiple root of the cubic equation at t = t*. The coefficients 
of this expansion depend on new scaling time r = {t — The leader 

term of the asymptotics and the corrections are defined by their asymptotics 

as r — oo uniquely. In particular, the leader term of asymptotics a (r) is 
a special solution of the Painleve-1 equation: 

S " _ 3 ^2 ^ ^ = 0. 
with the given asymptotics as r — > — oo: 

, . (5n-l) 1 1 

« (^) = 2^ "n^ ' > where ao = -y=, ai = — . 

n>0 V<3 ^4 

In the domain r > — oo this solution has poles on the real axis of r. Denote 
the least of them by tq. The asymptotics (|) is valid as (r — To)e~^^^ ^ 1. 

In the neighborhood of the bifurcation point (r = tq) the coefficients 
of the asymptotic expansion depend on one more fast time scale 6 = {t — 
To)e~^^^. Denote by 



n=l 

n 



where 6o= const. Then in the domain —e <^ 6'o <^ £ ^^^^ the formal 
asymptotic solution has the form: 

^{t. e) = (f/.+ & iOo) + e'l' E e^-'^l' w (^o)) exp (^) , 0. 
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The main term of asymptotics w (Oq) is the separatrix solution of the au- 
tonomous equation: 

tw' + U.(2\w |2+ + f/2 (^^*- ^) + I ^ |2 ^= 0, (6) 



namely: 







- Co) = T^r— cTp- 



In the domain — ^ 1 the asymptotic solution is defined by a sequence 
of two alternating asymptotics. Let us call them by "intermediate" and 
"separatrix" asymptotics. To obtain the intermediate asymptotics let us 
introduce one more slow variable: 

Tfc = ^fe-l£^/^A;=l,2,.... 

An asymptotic solution in the intermediate domain for not too large values 
k: 

k « 6-'/' 

has the form: 

^ CO 

^ n=0 

it 

The leader term of the asymptotics satisfies to the equation: 

Afc" + 3 Afc' = 

and can be expressed by the Weierstrass function with the parameter g2 = 0'- 

Ak=-2piTk;0,g3ik)), where g:,(k) = ^{gs{k - 1) + 7r/2), 93(0) = ^. 

(7) 

The constant 04 is the coefficient as (r — tq)^ in the Laurent expansion of 
a (r). 

The Weierstrass p-function has a real period Qk and has poles on the real 
axis at Tfc = and = Q^- The intermediate expansion with the leader 
term (|^) is valid in the domain between the poles as 

-e-'^'n > 1, e-^^'\Tk + Qk) > 1- 
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At the large values of k the intermediate asymptotics are constructed in 
the form 

/ °° 5n 

^ n=0 



2eJ' 



n=0 

The leader term of asymptotics satisfies to the equation: 

" „ , 
Ak +3 Ak = Afc. 



' j=l n=l j=l 

The main term of the asyptotics is: 

Ak (Tk) = -2p{Tk, Afc/2, gsik, e)), (8) 

where 

n °° 
g,{k,e) =9, (A;) + E^"/'° h (k). 

n=l 

The intermediate expansion with the leader term (^) is valid in the domain 
between the poles of the Weierstrass function as 

-e-i/^Tfc » 1, £-2/i^(Tfc + Qk) » 1. 

The separatrix expansions are valid in a small neighborhood of the Weier- 
strass function poles. Denote: 

= (Tfc + fifc - - E A; = 1, 2, ... . 

^ n=l 

When 

< 1 

the formal asymptotic solution of equation (|l]) has the form: 

^ = (f/*+ W {Ok) + e'" E W (^.)) exp (^) e - 0. 



The leader term of the asymptotics W {Ok) is a separatrix solution of the 



autonomous equation (H): 



-2 

W (Ok) = ■ 



'k-tU,)^' 
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The sequence of the alternating intermediate expansions and separatrix 
asymptotic ones is vahd as e~^^^{t^ — t) <^ 1. 

In the domain (t* — t)e~'^^'^ 3> 1 the asymptotic solution becomes two- 
phase. The amplitude of the stimulated oscillations in the solution of (|l|) 
oscillates fast. The form of the solution is: 



tlj= (^U iti,t,e)+eU (ti, t, e)+e^U (h, t, e)^ exp (^^^ , 
where ti is a new fast variable ti = S{t)/e + (f){t). The main term of the 



asymptotics U lies on the curve T(t): 

^\y\^-t\y\^-iy + y) = Eit), 
and satisfies to the Cauchy problem for the equation 





with an initial condition JJ \ti=o = uq, such, that 3(uo) = 0, ^{uq) 
minj^gr(t) (J^iu))- The function S{t) is a solution for the Cauchy problem 

^S' I = T, S\t=o = 0, 



m ^3y^ + {2E + t'^)y^ + 2ty + 1 

where T = const > 0. The function E{t) is the solution of the transcendental 
equation 



i / u*du = 71. 
Jr(t) 

The phase shift is defined by two constant. One of them is an "initial" 
condition 0(t*) = 00 and another one is a constant (pi in the equation in for 
0: 

Of.S 



Remark 1 In our work constants 0o ond 0i are not calculated. It means that 
the phase shift of the fast oscillating asymptotics remains to be undefined. Its 
evaluation within the framework of our approach requires an evaluation on 
an explicit form of the asymptotics of the higher corrections in the internal 
expansions at their matching with the fast oscillating asymptotics. 
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4 Construction of an algebraic asymptotics 



In this section the formal asymptotic solution for the equation (0) has been 
constructed in the form of series on integer powers of the small parameter e. 
We show that this expansion is valid in the interval {t — t^)e~^^^ ^ 1. 

The coefficients of this asymptotics are defined from the recurrent se- 
quence of algebraic equations. These equations can be obtained by substi- 
tuting of the series @) in the equation (H) and by collecting of the terms 
under the same powers of the small parameter. In particular, the algebraic 
equation @ for the leader term of the asymptotics can be obtained from the 
relation under 



Further we will use the root U2{t) as the main term JJ (t) of the asymp- 
totics d^). It means that we investigate the formal asymptotics (|D corre- 
sponding to the stable center which converges with the saddle at t*. 

The relation under give us the equations for real and imaginary parts 

1 

of the first correction term U= Ui + ivi 

[3 -t]Mi = 0, [U^ - t]vi = -dt U . (9) 

The expression [3 [/^ — t] vanishes at t = t^,. So it is easy to obtain that the 
even corrections are real, and the odd corrections are imaginary. 
The explicit form of the first correction is: 



1 i U 

U= -0 • (10) 

[3U'-t][U'-t] 

This expression has a singularity at t = t*. The order of the singularity 
increases with respect to the number of the corrections of the asymptotics as 
t — ^ t*. In the neighborhood of the point the coefficients of asymptotics 
(^) have the form: 

in the case of n = 0, 1: 

oo 

u{t) = it-Q--/'J2Ukit-Q'^^ (11) 

k=0 

^ 

where Uo= f/*, Ui= l/v3; 

in the case of m = 2n, m = 2n + 1, n G N: 

oo 

u{t) = {t-ur-'-^^/'j:uut-Q'/'. (12) 

k=0 
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The increasing of the singularity of the corrections in (^) makes this 
expansion unusable as t —>■ t*. The singularities of the coefficients define 
the domain of the suitability for the formal asymptotic solution as (t — 

Qe-^/" > 1. 



5 Expansions in bifurcation layer 

In this section the formal asymptotic solution for equation has been 
constructed in the domain |t — t*| ^ 1. Usually the asymptotics of this type 



is called internal ones |TH| 



5.1 Initial interval 

5.1.1 The first internal expansion — Painleve-1 interval 

In the neighborhood of the singularity we will use a new scaling variable 
T = {t — t^:)e~^^^. This new variable r is defined by the structure of external 
asymptotic solution singularities. 

We will construct the solution of the original equation (|]) in the form: 

U = U,+ e2/5a(r, e) + ie^'''[5{T, e), (13) 

The equations for a and (3 are 

a' + {Ul - U)I3 = -e^/'2U,a/3 - e^I\o? - r)/? - e^/^/J^^ 

/9' - W,o? + f/,r = e^l^{-o? - U,(3^ + ar) + e^'^a^^. (14) 

The asymptotics of the functions a and /5 as r ^ oo is known. It can be 
obtained by substituting (|ll]) and (0) in formula (|[) and by decomposing 
of this expression in the terms of the scaling variable r: 

1 oo OO ^ 

V 3 n=l n=0 

4/5 rr„ ^(3-5n)/2 



e 

oo „ , 1 oo 



E t/2 r(3-^")/2 + . . . . (15) 



n=0 

1 ^ 2n+l n I K >/o o/K 2n+l 



n=l n=0 



It is convenient to construct the solution of the equations (|14D in the form 
of formal series on powers of the small parameter e: 



E 3 (r), /3(r, 5) = E ^ (^)- (17) 



n=0 n=0 
11 



The leader correction terms are defined by the system of ordinary differential 
equations 











a'+{U',-Q P=0, (18) 
^'-3f/, a^ + rf/, = 0. (19) 

Let us differentiate the first equation and take into account that U^,{Ul—t^) = 
1 then we obtain the Painleve-1 equation for a 

a" -3a^ + T = 0. (20) 

Matching the asymptotics (p!3|), ( p!7D and asymptotics (H) as t — + we 
obtain the asymptotic structure of a (r) as r — > +oo 

" (^) = a„r"^^, where ao = 4^, = ttt- (21) 
n>0 v<3 ^4 

The solution of the Painleve-1 equation with the asymptotics of ( [21| ) was 



investigated in [17, 18 



The coefficients of the formal expansion ( P^D for the function (3 can be 
obtained from the first equation of the system (0). In particular: 

a' 1 1 / 1 n \ 

2 1/2 nO iO nl o\ 

/3= ^ a'+ /5 +2[/, a/3 +2?7, a/3 -r /3 j . 
The equation for the first correction term a (r) has the form 
a " + 6f/,(f/2 - Q aa +2f/,(f/2 - t,) + {U^ - U) a' - {Ul - U) a t = 

-2U,a'/3 -2U, a^ '. (22) 

The asymptotics of a solution for this equation is uniquely defined by match- 
ing of the asymptotics of the external expansion (^ and internal expansion 

(0), m-- 

a=U2 r- Ue t'''/^ + 0{r~^), r ^ +oo (23) 

A simple calculation gives us the domain of a suitability of that asymptotic 
expansion. The representation (|l^ of the solution is valid as r <^ e~'^^^. 
This result is similar to [|5[. 
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It is known that the solution of the Painleve-1 equation with the asymp- 
totics (|2T|) has poles at points r^. Let us denote the first pole by tq and 
call it the point of the bifurcation. 

The solution of the Painleve-1 equation is represented by converging series 
in the neighborhood of the point Tq | T8| : 



2 T °° 

" (^) = -7 ^^ - T7;(^ - ^0)' + ^^(^ - ^0)' + J2^kir- ro)^ 

(r - Toy 10 

the constants tq and 04 are parameters of the solution. 

The first correction term a is defined by the equation ( ]2^ ) and asymp- 



totics (]23|). In the neighborhood of the bifurcation point tq as r > tq the 



function a is represented in the form: 

a=ac a+ai(r)+ 6^a2(r). (24) 

Where (r) is the partial solution of the nonhomogeneous equation for the 
function a (r) with the asymptotics: 

11 ^1 ill 

ttc (^) ="-4 (r - To)"^ + ^ ftfc (r - To)'', a_i= 04=0. 

A:=0 ^* 

It is convenient to use the partial solution with the coefficients as (r — tq)~^ 
and (r — Tq)"^ equal to zero. 

The constants a+ and 6"^ in the formula (^) can be uniquely defined by 
asymptotics of a as r ^ 00. The functions ai(r) and 0:2(7") are linear inde- 
pendent solutions of the homogeneous equation for a (r). These functions 
are uniquely defined by their asymptotics as r — > tq: 

ai(r) = (r - Tq)"^ + ai(r - tq) + a2(r - tq)^ + a^ir - Tof + 0((r - Tq)^), 
where a„ = const; 

«2(r) = (r-ro)^ + 0((r-ro)«). 
The next correction terms are represented by the following formula: 

71,, "/ \ "4- / \ "4- / \ 

a (r) =a^ {t)+ a+ai(r)+ 6+02(1")- 
The partial solution Sc (r) of the equation has the asymptotics as r — > tq: 

Oc (r) =a-2n (^ - ^0) + • • • , 
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where dots mean the terms of the next orders on (r — tq). 

It is convenient for us to denote the partial solution which does not contain 
terms of the order of (r — Tq)"^ and (r — Tq)^ in the asymptotics as r ^ Tq 
by Sc (r). 

By analysis of written asymptotics we obtain that the formal asymptotic 
solution (|l^) is valid as (r — ro)£^^/^ ^ 1. 

5.1.2 The second internal expansion, the neighborhood of the bi- 
furcation point 

In this section we construct the formal asymptotic solution for the equation 
(0), which is valid in the neighborhood of the bifurcation point |r — tqI <^ 1. 

The domain where the first asymptotic expansion is valid determines the 
second internal scale. We will use a new scaling variable 9 = (t — ro)e~^/^ in 
the small neighborhood of the point tq. 

We change the unknown function in the equation and pass to the 
function w = U — U^. 

Substituting the expressions for w and 6 in the original problem we obtain 
the equation for w{9,e) 

iw' + 2f/> + 2U,\w\'^ + \w\'^w + Ulw* + U.^w'^ - 

-e^'^TQiU, + w)- e0{U, + w)=0. (25) 

The asymptotics of the function w{6,6) as 6 ^ +oo can be obtained 
by expansing in series of the function U = + €'^^^a{T,e) + ie^^^P{T,e) on 
variable 9: 

£3/5 + ...)+ (^^c' + ...) + .... (26) 

In this formula dots mean the terms of the lesser order (in brackets with 
respect to 9 and out of them with respect to e). 

In this section the formal asymptotic solution for the equation (^) with 
the given asymptotics (^) as e ^ has been constructed. We indicate the 
domain where this formal solution for the equation (^) is valid. 

The first step of the formal solution construction for the equation (|25| ) 
is a choice of an asymptotic sequence. This choice for the external and the 
first internal expansions was defined by the equation. In this case the asymp- 
totics as 6* — >• +00 gives us the asymptotic sequence on e^^^. But the equation 
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shows, that the first order correction should satisfy to the homogeneous equa- 
tion. In this case we can satisfy to the asymptotic condition as 6 —>■ +00 by 
constructing of the formal solution in the form: 



00 



(27) 



n=l 



where 



e + e^l^ On 00+... e 



0. 



The correction terms are defined by the asymptotics of the formal solution 

as ^ 00: 

n , 



It is convenient to pass from the equation 
to the variable 



to the equation with respect 



iw' + {2Ul - t,)w + 2U^\w\^ + \w\^w + Ulw* + U^w^ = 



(f/* + w) 



{2i 



To construct the formal solution w{dQ,e) we substitute the formal series 
(|27|) in the equation (^) and gather the terms of the same order on e. Then 
we obtain the recurrent sequence of the equations for the leader term and 
the corrections of the asymptotics. The equation for the leader term is 



iw' + UA2\w\'^+w'^]+ Ul [w*- w ] + \ w l"^ w= 



(29) 



The asymptotics of w as ^0 ^ 00 can be obtained from the asympotics 
of the function w{9o,e) 



1 



(30) 



To construct the leader term w (9q) we take the conservation law for the 
equation (p9|): 



2^ ' 2 ' ' . 



1, , 
+ 2'^' 



H. 



The constant H for our solution is defined by asymptotics (|30D: H = 0. Let 
us solve this formula for the function w* and substitute the representation 
over w in the equation (^). Then we obtain 



iw' — 2U^w\ —U^w = 0. 



(31) 
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This equation can be obviously solved 



a (««) - 7^^. (32) 



the integration constant is defined by the asymptotics (pOf) . The leader term 
of the asymptotics ( pTl) is the separatrix solution of the equation (pQ]). 
The next corrections w, m = 1,2 . . . satisfy to the equations 



i w + 



2f/, w* + 2f/, w -Ul + 2 WW* 



/ n \ 2 



w + {U,+ w\ w*=F. (33) 



In particular, 



F=To{U^+w), F=eo{U^+w) F=-0o{U^+w), F= 



F= - ^ (f/*+ ^) - 2f/,| w p - f/, 

The asymptotics w (6'o) as — oo can be easy obtained from 

To construct the correction terms we need the solutions of the homoge- 
neous equation (|3^) . One of them can be obtained by differentiation with 

respect to variable of the function w 

The formula for the second independent solution is complicated and we 
do not write it. The second solution grows as |^| — > oo. The asymptotics as 
6' — > oo of this solution has the form 

w2{eo) = el + i^^el + .... (35) 

These independent solutions are such that their Wronskian is 

7i 

W = wiwl - W2WI = - — . 

u: 

By using the expressions for wi and W2 we can write the next corrections 
for the asymptotic solution with the given asymptotics as +00 in an 

explicit form. The asymptotics as 6'o —>■ — cxd of these corrections can be 
evaluated. 

It is convenient for us to represent the higher corrections of the asymp- 
totics (|27|) in the form 

n ri- ^ n 

w (9o) =w, {9o)+ Xo~Wii9o)+ Yo'W2i9o), (36) 
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where, Wf. {Oq) is the partial solution for the n-th correction term and this 
solution does not contain terms of the order of 6q and Oq^ in the asymptotics 

n n 

as Oq = oo. The constants Xo~ and Yo~ can be evaluated by 



-reSr 



— resr 



- jJFw,+ Fw,)- 



The segment of the Laurent series for the function w^. (9q) has the form 



1 



5 " 5 30U, 
The coefficients as Wi and W2 are 



0% 



oo. 



1 



0, Fo" 



0. 



The function w (9q) a,s 9q ^ —oo also can be represented in the form of 
sum of the partial solution of the nonhomogeneous equation and solutions of 
the homogeneous equation 



-Wr 



(37) 



The partial solution of the nonhomogeneous equation can be represented in 
the form as 



2 

Wr 



6 



.U,9l 59o .3 

I H i- 

2 4 



15^7, 



16f/2 2 



15 



3f/*' 



0(V 



The Laurent series of this solution does not contain the terms of the order 
of 9q and 9q^. The coefficients as wi and W2 are 



2 



0, 



2 



TT 

2 



The next correction terms of the asymptotics (|27|) also can be obtained 

by this way. The condition e w / (e^^^ -C 1 gives us the domain where the 
expansion ( p7|) of the asymptotic solution is valid 



-1/5 



< 6^0 < £ 



'1/10 
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Remark 2 It is easy to see that the asymptotics of as 60 —>■ oo and as 
6q —00 are different. The difference is the term of the order of eO"^. It 
means that additional terms appear via the passage of the neighborhood of the 

2_ 

bifurcation point. The leader term of these additional terms is Yq = 6119 /2. 

This effect leads to the appearance of the solution Yq W2 of the homogeneous 
equation in the asymptotics of the second corrections as Oq —00. 

5.1.3 Intermediate expansion, initial interval 

In this section we investigate the behaviour of the solution after of 
the narrow neighborhood of the bifurcation point Tq. 

In this domain it is convenient to introduce a new scaling variable Ti = 
9q£^^^ by analysis of the second internal expansion. 

The solution of the equation we seek in the form: 

U = U, + £2/6a(Ti, e) + ^£^/^5(Ti, £), (38) 

where A{Ti,6) and B{Ti,e) are real functions. 

The equation (j^) can be written in terms of these new variables in the 
form: 

A' - 2UlB = -e^l^2U^AB - e^'^A^B - eB^ + 

00 

(^4/5^^ + ^5/6^^ _ ^4/5 ^ ^(2n-l)/5 

n=l 

B' + W,A^ = -e^/'\A^ - B^U,) - e^/''AB^ + {e^'^T^ + ^^/i^o - 

00 

^2/15 ^^(2n-l)/5 +^1/3^)^^ ^ Q. (39) 

The asymptotics of [/ as ^ — oo rewriting in terms of the variable Ti 
gives the asymptotics of A and 5 as Ti — —0: 

A = (^+ Fo T,^ + ...)+ e^/^^ {^T^ + ...) + 

.^/^(^T3 + . ..) + ..., 
i?=(^+roTf + ...)+.^A^(^T, + ...)+.... (40) 

The dots in the brackets mean the terms of the less order with respect to Ti 
and out of the brackets - with respect to e. 
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The system ( |39| ) and the asymptotics ( ^0]) contain terms of the different 
orders with respect to e. It is convenient to represent the formal asymptotic 
solution with respect to e in the form: 



oo 



n=0 

oo 
n=0 

By substituting the representation for A and B in the equation we obtain 

n n 

the recurrent sequence of the equations for A (Ti), B [Ti). It is convenient 
for us to investigate the second order differential equation for the functions 

n 

A (Ti) instead of the system of differential equations of the first order for 

n n 

A (Ti), B (Ti). The equation for the leader term is 

A" + 3A^ = 0. (42) 
The equations for the next corrections of the asymptotics (^Tj) are 

in particular, 



m 0mm 

A" + Q AA=F, 



F=0, m = 1, 2, 3, 6, 7; F= -2f/,Vo; F= -2U!T,; 

8 4„ 9 4 5 

F= -3 aHTi); F= -6 A (Ti) A (T,); 

9 - 

10 1 /( AT 3 A\ o/ o„ (A') 



2U, 2U^J *\ 4[/3 

The asymptotics as Ti — > —0 of the coefficients of the formal series for 
A[Ti,e) can be obtained from the formulas P0| ) 

0—22, m 

-4= Y^Tf + . . . ; yl= 0, m = 1, 2, 3, 6, 7; 

5 
l(Ti) = 0(Tf); 1(Ti) = 0(T7); A (Ti) = 0(Tf 

The dots in the formulas mean the terms of the less order with respect to 
Ti. In particular, these formulas give us the domain where the segment of 
the length n = 10 of the formal series (^if ) as Ti ^ — is valid 

£i^^«l, orT,».Va 
A{T,) 
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Let us construct the formal asymptotic solution in the form (^T]). The 

seeking solution A (Ti) for the leader term equation (|^) has the asymptotics 
as Ti —>■ —0 in the form 

A=-^ + a,{l)T^ + 0{Tn, 

2_ 

where 04(1) =Yo ■ Using this asymptotics it is easy to see that the function 



A is expressed by the Weierstrass function 

l(Ti) = -2p(Ti; 0,^73(1)), where g,{l) = a^{l)/5Q. 

n 

In a general case the solution A (Ti) can be represented in the neighbor- 
hood of the point Ti = in the form 

A (Ti) =M{Ti)+ Xo~Ai{T,)+ WA^rn. 

n 

The function A§{Ti) is the partial solution of the nonhomogeneous equation 
for the n-t\i correction. This solution does not contain terms of the order 
of T{'^ and Tf. The functions AiiTi) and y42(Ti) are linear independent 
solutions of the linear homogeneous equation in variations. These solutions 
are uniquely defined by their asymptotics as Ti —0 

A(Ti)=Tf3 + ..., A2{T{)=T^ + ..., 

the expansion of the function y4i(Ti) does not contain terms of the order of 
T^. Both of these solutions can be represented via the Weierstrass function. 
In the case = we obtain (PU[) 

Ai(Ti) = ^p'(Ti, 0,(73), 



A2(Ti) = -14c',3p(Ti, 0,(73) = r^Tip'(Ti, 0,(73) + ;J^p(Ti, 0,(73). 



The Wronskian of these solutions is W = — |. 

n n 



The coefficients Xq Vq are defined by matching of the asymptotics (HTl) 
as Ti —0 and the asymptotics ( pT] ) as ^ 00 



n 



Xq =Xq , n = 6m + 23; 

n m 

?/o~ =Fo~, n = 6m-12. 
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u 

The function A has poles. One of them is Ti = 0. The second pole of 
this function is defined by an elliptic integral 



dy 



The formal asymptotics (0) loses the asymptotic behaviour in the neigh- 
borhoods of poles. The asymptotics as Ti —0 is known. These explicit 
formulas give us the domain where the asymptotic series (^1]) is valid as 
-Ti > £^6. 

To construct the asymptotics of the coefficients of the formal solution as 
Ti —Qi + we solve the equations for the correction terms and write their 
expansions as Ti —Qi + 0. 

n 

It is convenient to represent the correction A as Ti ^ —Qi + in the 
form 

A (Ti) =Ac (Ti)+ x^+A^{T^ + yi+A2(Ti + n,). (43) 

n 

The function Ac {Ti) is the solution of the nonhomogeneous equation for the 
n-th correction. The asymptotics as Ti —Qi + of this function does not 
contain terms of the order of (Ti + fii)^'^ and (Ti + fii)^. 

The constants in the formula for the n-th correction are evaluated by 
formulas 

n 

"1 r-^i+r dz F {z)Ai{z) 

_r J-r 



n , 
Xi 



--Xq + Ci + reSr=o 
where the constant Ci is defined by 



W 

n 

1 /--f^i+r- dz F {z)A2{z) 

r J-r 



W 



A2{Ti + n,) = CiAi{Ti) + A2{Ti). 

To obtain the domain of the suitability of the formal asymptotics as 
Ti —Qi + it is necessary to know the type of the singularity of the 

4 

correction terms. The first nonzero correction term is A (Ti). The expression 

4 

(|43|) for n = 4 contains the partial solution Ai(Ti). This function is smooth 
and has the identical asymptotics as Ti —* —Qi + and as Ti — >■ —0. But 

4 

the asymptotics of the A SiS Ti —Qi + contains terms of the order of 
(Ti + fii)^ with the coefficients Xi' ^ 0. Namely, 

4 , 

= 0, 
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where ({Ti, g2, gs) is the Weierstrass ^-function. 

The values of the coefficients Vi^ and Sf*" for n > also are evaluated. To 
obtain the domain where the asymptotics (^T]) is valid it is necessary to note 



It is easy to see that the condition e~'^^^^{^li + Ti) ^ 1 is the condition of 
the suitability of the asymptotics (|4ll) as Ti — > —Qi + 0. 

Thus, we have constructed the formal asymptotic solution of the origi- 
nal equation in the form ( PD in the domain where the following conditions 

-e-^/^Ti > 1 and e-^/^^i^i + Ti) > 1 are vahd. 



5.2 The bifurcation layer in the case of bounded k 

5.2.1 The neighborhoods of the poles Tk of the Weierstrass func- 
tion 

In this section we construct the solution of the original problem in the neigh- 
borhoods of the poles of the solutions of (|42|) and show the change of the 
solution via the passage from one pole to another one. 

In the neighborhood of the poles we will use a new scaling variable 

1 oo 

e, = {T, + n,--Y: E ^ 0. (44) 

^ n=l 

The original equation (§) looks like 

/ k -| oo k \ \ 

( E + I E E 5/ + ^""^^ p = i> ^ 0, (45) 

^i=l ^ n=\ i=l ^ j 

where Tq is the bifurcation point and fi^ are poles of the Weierstrass func- 
tion which describe the leader term of the asymptotics in the intermediate 
domains. 

The solution is seeking in the form of formal series 

°° n 

= ly (4) + 5^/5 (^fc)' ^ 0- (46) 

n=l 
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The coefficients of the formal expansion (^Uf ) are satisfied to the recurrent 



system of equations. In particular, the function W (Ok) is satisfied to the 
equation (29). The equations for the higher corrections of the asymptotics 
are 



i W 



2U, W* + 2U, W -Ul + 2 WW' 



\ 2 n 



W + [U,+ W] W* =F ■ (47) 



n m 3 \ 

Let us become clear how the nonlinearity F depends on W and on X]^ . 

n 

The maximum order of nonlinear terms is cubic. The function F contains 



terms of the type 



and 



j I m 

WWW, n = j + l + m + AQ, j,l,m>l; 



{U^+W)WW, n=j + l + 23,j,l>l] 



linear terms with constant coefficients with dependness on x]^ 

k 

W^Xi, n = m + j + 25, /, m,j > 1; 
1=1 

linear terms with the coefficients depending on 6k 

m 

Ok W, n = m + 55] 
and with coefficients only depending on x^ 

^ 

^x/(f/,+ Vr), n = m + 2; 

m 

and linear terms with depending on 6'o- 

m I 

OoW, n = 12m + / + 18; 
and the terms depending on the leader term of the asymptotics 

m 

eo{U*+W), n = 12m-5. 

For example, we write a few of the right hand-side terms for the lower 
corrections: 

1 2 JL n 

F= To{U,+ W); F= (f/*+ W) J2 F=0, n = 3, 4, 5, 6; 
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"1 k 

F= {u.+ wm+jY^'^n 

It is easy to see that the right-hand side of the equations for a few of the 
lower corrections are proportional to leader term of the asymptotics. The 
right-hand side for the 7-th correction term depends on the leader term by 
more complicated way. There is the variable coefficient 6^ in the equation. 

The solutions of the equations for the corrections terms of the asymptotics 
can be represented in the form 

71 71 71 TL 

W {9k) ^Wc (Ok) + X^wi{9k) + Yk^W2{9k). (48) 

n 

Where the function Wc {9k) is the solution of the nonhomogeneous equation 
for the n-th correction term. This solution does not contain the terms of the 
order of 9f. and 9^^ in its asymptotics as — cxd. The linear independent 
solutions of the homogeneous equation Wi{9k) and W2{9k) are defined in the 

section 5.1.2. The constants Xk and Yk will be determined below. 



The asymptotics of the function W {9k) as — > cxo is given. It can be 
obtained by rewriting of the formal solution + s'^^^A{Tk, s) + is^^^B{Tk, e) 
as Tfe — > —ilk + in the terms of the scaling variable 9k. This asymptotics 
uniquely determines the leader term of the formal solution 

Thus, the dependence on the number k in the formula for the leader term of 
the formal series is only contained in the argument of the variable 9k- 

n n 

To evaluate the constants XiT and Yj^ in the formula for the n-th correc- 
tions we match its asymptotics as — > oo with the asymptotics obtained by 
rewriting the formal series + e'^^^A{Tk, e) + is^/^B{Tk, e) asTk^ -Qk + 
in terms of the scaling variable 9k- We obtain 

= 0. (50) 

The reason of this effect is that all terms connected with the solutions of 
the type of a;fc+y4i(T/j) have been matched by the shift of the independent 
variable 9k- 

To obtain the formula which connects the constants Yk'^ and Vk'^ it 
is necessary to note that solutions A2{Tk) is matching with the solutions 
W2{9k)- Thus the following formula is valid 

71 TfL 71 

Yk^=yk^, + = 0, n = 0, 1, . . . , 6; (51) 
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7 

The value of the coefficient Yk is determined by matching with the asymp- 



totics of A (Tk) 

Yft = 56(73 (A:). (52) 
Let us represent the expansions of two corrections of the (|^) as 6'^ ^ oo: 

W {9k) =Wc (Ok), (53) 

where 

W. {9k) = -^9, + ^^9k - ^ + ^^,9, + 0{9, ) 

These formulas explicitly contain the bifurcation point tq. The terms of the 
asymptotics begin only to depend on the poles Qk, k = 2, ... of the Weier- 

2 

strass p-function with the second correction term W {9i). This correction as 
9k —>■ oo is represented in the form 

W {9k) =Wc {9k) (54) 

where 



2^ 



W. {9k) = -^9, + ^^9k - ^ + ^^9, +0{9, , 

k 

Pk = ^Vtj, 9^^oo. 
i=i 

The formulas (PI), (^6D, (^)- ( ^2]) determine the structure of the formal 
solution in the neighborhood of the point 9k = +oo. It is easy to see that for 
Vfc the formal solution (^) is valid as 9kS^^^ <^ 1. 

Now we know the form of the solution as ^ oo and two linear in- 
dependent solutions of the equation in variations and we are able to write 
the representation of the correction terms in the formal solution (^61) as 
9k -oo. 

n Tl 71, Th 

W {9k) =Wc {9k) + Xk wi{9k) + Yk W2{9k). (55) 



Where the functions Wc {9k) have the similar form as in the formula (^ 

n n 

The constants Xk and Yk are 



Xk =Ak , Afc = -res^ 



1 /"^ n ^ ^ d9k 

-JjFw;+F*W2)^ 



(56) 



25 



Yk =Yk+ Ak , Ak = -resr 



1 Z""" / " * N d9k 



(57) 



It is important to note that 

n 

n- = 0, as n = l,...,6; (58) 
Yk =n+ + |. (59) 

Remark 3 The formula (|5£) in the case of the second order differential 
equation was obtained in the work j^. 



The domain where the formal series ( ^61) as 9^ — oo is vahd is deter- 
mined as well as — *• cxD. We obtain that for any value of k this domain is 
determined by —OkS^^^ -C 1 as well as at large values —Ok- 

5.2.2 The intermediate expansion, the internal interval — interval 
of the Weierstrass p-function 

In this section we study the formal asymptotic solution in the domain be- 
tween k-th and k + 1-th poles of the Weierstrass function. The solution in 
this domain has the asymptotic structure which is similar to the asymptotic 
structure in the initial intermediate interval. The difference is: the asymp- 
totics on the right side of the investigated interval is matching with the 
asymptotics (^) where the asymptotic sequence is {e"'^^^}. But the formal 
asymptotic solution has been constructed in the intermediate domain on the 
initial interval and has been matched with the asymptotic sequence {e"/^}. 

In this intermediate domain we use the slow variable Tk = OkS^f^ and seek 
the solution of the equation (|^) in the form 

U = + e^/'Ak{n,e) + te^/'Bj,{Tk,e). (60) 
The equation (H) becomes 

A', - 2UlBj, = -s'/'2U,AkBk - s^'^AIB^ - eBl + 

oo 

■^4/5^Q_^4/5^^(2n-l)/5 ^ 



n=l 

k oo k \ 

e''\n - E - E E 

j=l n=l j=l ^ 

B', + W.Al = -e'/'{Al - BlU,) - e"'A,Bl + 



^2/15^^ _^ ^2/15 ^(2n-l)/5 



n=l 
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The asymptotics p6| ) as — >■ — cx) of the function U has been written 
in terms of the new variable gives us the asymptotics as — > —0 of the 
functions A and B 

A, = (^+ Y +...)+ e'"' {^-^n +...) + 

= (^+ F^_,T3 + ...)+ s^A5(lflT, + ...)+.... (62) 

The dots in these formulas mean the terms of the small order with respect 
to Tfc in the brackets and with respect to e out of them. 

k 

The system of the equations (^) and the asymptotics (|62|) contain terms 
of the different orders with respect on e so we seek the formal asymptotic 
solution in the form 

oo 
n=0 

oo 

B,iT,,e) = ^£-/30 (T,), e -> 0. (63) 

n=0 

Using the anzats for the functions and -B^ we obtain the recurrent 

n n 

sequence of the equations for Ak {Tk)-, Bk {Tk)- As it was done above it 
is convenient for us to investigate the second order differential equation for 

Ak [Tk) instead of the system of the first order for Ak [Tk), Bk [Tk)- The 
equation for the leader term of asymptotics ( |63l) is 

Ak" + ^Ak^ = 0. (64) 
The equations for the next correction terms are 

m m m 

Ak +QAkAk=F . 
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The expressions for the right hand terms F depend on n. Here we represent 

n 

the dependness of the F on corrections terms with a small numbers and on 

J + n 

the coefficients . The right hand functions F contain terms of the type 

m / 

with quadratic dependence on corrections Ak and Aa; as n = m + / + 10; 

m j I 

with cubic dependence on corrections Ak, Ak and Afe as n = m + / + j + 20; 

m j + 

with linear dependence on corrections Ak with coefficients depending on 
as n = m + j + 15; 

m 

with linear dependence on correction Ak with the linear coefficient with re- 
spect to (Tfc + Z]j=i ^fe) as n = m + 15; 

with linear dependence only on as n = j + 5. 
In particular, 

F=0, m= 1, 2, 3, 6, 7; F= -2U^to; F= -2U^Tk + Pk] 
F= -3 Ak\ny, F= -6 Ak in) a (Tfc); 

The asymptotics as ^ —0 of the coefficients of the formal series can 
be obtained from the formulas (|62|) 



0—22. m 

Ak= -7^+ YkT^ + ...; Ak= 0, m = 1, 2, 3, 6, 7; 

Am)=0(T«); 1*=0(T«); ^1= 0(7^7*). 

In these formulas the dots mean the terms of the less order with respect to 
Tfc. In particular, from these formulas we obtain the domain of suitability 
for the segment of the series (|^) as Tk —0 

^^-^ « 1, or -.-V6T,»1. 



Ak (Tk) 

Let us construct the formal solution in the form (|63|) . The seeking solution 
for the leader term of the asymptotics equation should have the asymptotic 
structure of the type as — >• — 

Ak=-^+Y~,.,T^ + 0{T^'). 
28 



Using this formula it is easy to see that the leader term Ak of the asymptotics 
is expressed by the Weierstrass p-function 

Ak {Tk) = -2p{n; 0,g3{k)), where g^ik) =Yk-i/5Q. 

n 

In a general case the solution A (T^) can be represented in the neighbor- 
hood of the point = in the form 

Ak {Tk) =M{Tk)+ XkA,{Tk)+ VkA^in). 

n 

The function Ak'^{Tk) is the partial solution of the equation for the ra-th 
correction term. This solution does not contain terms of the order of Tj^^ 
and in the asymptotics as — > 0. 

n n 

The coefficients Xk' and Vk' are defined by matching of the asymptotics 
(|63|) as Tk —0 and the asymptotics ( ^61) as 6*^ — oo 

n _ m _ 

Xk =Xk , n = 6m + 23; 

n m 

Vk =Yk , n = 6m - 12. 



The period of the function Ak {Tk) is defined by the elliptic integral 

^4^3 -g^(k) V 4 ; 

The formal series solution is not valid in the neighborhood of the 
poles. The asymptotics as ^ of the coefficients of the formal solution is 
known. These explicit formulas give us the domain where the formal solution 
(|6|) is vahd. This domain is -e'^/^Tk > 1. 

To construct the asymptotics of the coefficients of the formal series solu- 
tion in the small neighborhood of the pole = —Qk we solve the equations 
for the corrections and rewrite their expansions as —>■ —Qk + 

n 

It is convenient to represent the correction term Ak as — * —^k + in 
the form 

Ak (Tk) =Ac {Tk)+ Xfc+i+Ai(Tfc + nk)+ Vk+^A^iTk + fifc). (65) 

n 

Where the function Ac {Tk) is the solution of the nonhomogeneous equation 
for the n-th correction and its expansion does not contain terms of the order 
of {Tk + Vtk)~^ and {Tk + Vtk)'^ as Tk -Qk- 
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The constants in the formula for the n-th corrections are evaluated by 



n , n n n 

Vk+i =yk+ 5^ = res,=o 



1 f-^k+r dz F (^)^l(^) 



n_|_ n_ X X 

^k+i + t^fc y k+1 + Sk , 6k = reSr=o 

where the constant Ck is defined by 



r J-r W 

n 

1 f-^k+r dz F (-2)^2 (-2) 



W 



A2{Tk + Qk) = CkAiiTk) + A2{Tk). 

To determine the domain of the suitability for the formal series (|63D as 
Tk — i> —^k + it is important to know the structure and orders of the 
singularities of the correction terms. The first nonzero correction term is 

4 4 

A (Tk). The representation ( ]65|) for n = 4 contains the function Ak+i(Tk) 
as the partial solution of the correction equation. This solution has the 
identical structure of the asymptotics as Ti ^ fii + and as Tk —0. But 
the coefficient Xi' 7^ as Ti ^ fii + 0. Namely, 

k+i^ = 0, 
^k+i' = ^^C(^fc/2;0,^3), 

where ({Ti, g2, Qs) is the Weierstrass ^-function. 

The values of the coefficients Z/f^ and x^ for n > can be also evaluated. 
To obtain the domain of the suitability of the formal series (p3|) we note 

5^/15 A (Tk) e'/'' 

— -const <^ 1. 







A{T> 



k 



(Tk + ^k) 



It is easy to see that the condition e~^/^^(Tfc + Vlk) ^ 1 describes the domain 
where this solution (|63|) is valid as Tk —^k- 

Thus in this section we have constructed the formal series solution in the 
form ( |63D and have shown that this representation of the solution is valid as 
-e-^/^Tk > 1 and e'^'^^Tk + fi^) > 1. 



5.2.3 The discrete dynamical system 

The coefficients which determine the intermediate and the second internal 
solution are changed by passage from the A; to fc + 1. Let us investigate the 
discrete dynamical system for the coefficients Xkt Yk and 6*^, x^", 2/^ . 
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The connection formulas for these coefficients are 

y k+1 =Y k + 

where n = 6m — 12, 

m _ n _|_ m y 

Yk+1 + Afc+i- 

m 

Where Yk= as n = 0, . . . , 6. In the case of m = 7 we determine the 
parameter (73 (A;) of the Weierstrass p-function for the leader term of the 
intermediate expansion 

g3{k) = ^ia, + ^ik-l)). 

The connection formula is defined by explicit way without investigating of 
the next corrections because at first 



6k' = 0. 

and secondly 

Remark 4 This connection formula for g^lk) in the case of second order 
ordinary differential equation was found in the work /0/ without investigation 
of high correction terms in internal domains and discrete dynamical system. 

It is easy to see that this system for Y k+i Uk+i is not closed. It 

is connected with the fact that the terms A ^+1 and S k+i contain the terms 
I _ j 

Xk and Xfc . So to obtain the closed system for the coefficients two above 

I _ j 

relations have to supplement by equations for Xk and . As result we 
obtain 

^k = Ck Vk + 6 k-i + A k~i- 

This discrete dynamical system allows us to evaluate the phase variable 
in the leader term of the second internal expansion 

9k+i =9k + e-^'^ (fife+i + ^ f; £"/3o 5+^^ \ ^ ^ 0. 

^ ^ n=l ' 

Thus, to obtain the phase variable to within o(l) in the neighborhood of 

the [k, + l)-th pole it is necessary to obtain the term x^+f" as n < 5 and Y k 
as m < 3. 
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Let us consider these formulas as n < 5 and m < 3 

Afe" = 0, m 7^ 7; 

n 

= 0, < n< 5; 

I 

We also have that Afc = 0, as / < 1. 

As result the formulas for Xj^ can be written in the form 

= = 0, n = 1, 2, 3, 

Z/fe = 0, Xfe =(5/f , 

5.2.4 The domain of the suitability of the intermediate expansion 

The domain where the intermediate asymptotic expansion is valid for big 
values of k is defined by inequality 

^ J^- 

Ak 

5 5 

The function Ak contains the term Vi^AiiTk) which is the biggest with 
respect to k 

This formula leads us to the inequality 

£^/6A;^/6 1, A;«£-^/^. 
It is easy to determine the order of as /c — > oo 

k 

J=l 

5.3 The intermediate expansion for large k 

In this section to take into account the secular terms as large values of k we 
modify the intermediate expansion. 
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5.3.1 The domain for moderate large values of k 

As was shown above for the large values of k the correction terms of the 
asymptotics contain the secular terms with respect to k. These singularities 
have the different structure. At first, the singularities are connected with 

n 

the terms Uk as the solutions of the homogeneous equations. Secondly, 
the singular terms Pk and Y,^=i are contained in right-hand sides of 
the equations (|6lD. These singularities lead to lose of the suitability of the 
equations for the correction terms equations. 

To suppress these singularities we allow the dependness of the parameter 
Qs on the small parameter 

n °° 
gsik,e) =9s ik) + Y.e^/''' 9^ (k). 

n=l 

Where 

93{k + l) =93{k) + —{Ak + 6i), m = n + 5. 
The equations for become 

Xk^ = Ck 5k^ + , n = / + 29. 

We can suppress the first kind of the singularities in the intermediate expan- 
sion by this way. To suppress the singularities of the second type connected 
with large values of Pk and sum of the coefficients Xj^ we include the following 
term in equation for the leader term of the intermediate expansion 




Then the leader term equation has the form 

Ak +3 Ak = Afc, 

and the solution is 

Ak (Tk) = -2p{Tk,\k/2,93{k,e)). 

Including of the term Xkis) in equation for the leader term of the asymp- 
totics leads us to change of the asymptotic sequence. Namely, 

00 5n 

A{Tk,e) = Y.e^/^ A{Tk). 

n=0 
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The terms of the order of e"^/^^ with m 51,1 E Z are contained in the 
leader term of the asymptotics and corrections by introducing of the Xk{£)- 
To obtain the intermediate expansion with respect to should expand 

this representation with respect to e. 

The solutions of the corrections equations are defined by their asymptotics 
as Tk -0: 

5n 5n 571 

Ak {Tk) =Ak%Tk)+ xfA.iTk) 

The main term of the asymptotics with respect to Tk is the term of the 

order of T^. The function Ak is the partial solution of the nonhomogeneous 
equation for the correction term and does not contain the terms of the order 
of T^ and T^7^ in the asymptotics as —0. The function Ai is the 

solution of variation equation and is defined by the leader term 

A^ = -^p{n,Xk/2,gs{k,s)). 
aik 

In the neighborhood of the point T^ — — fife + the corrections can be 
represented by formulas 

Ak (Tk) =MiTk)+ H+iA^iTk) + llA^iTk). 
The function A2 is the second solution of the homogeneous variation equation 

A2 = ^^(^fc' W2, gz), g3 = gsik, e). 

The constructed solution is valid between the poles of the Weierstrass 
p-function p{Tk, Afe/2, g3{k, e)) in the similar domains as in the case of inter- 
mediate expansion for finite values of k. 

5.3.2 The second internal expansion for moderate large k 

To construct the second internal expansion in this case we should make the 

n 

changes in the formula for evaluating of YjZ- 

= 56 (A;), n = 6m - 12. 

5.3.3 The large values of k 

The large values of k lead to grow of the term Afe(e) in the right hand side 
of the equation for leader term of the asymptotics. The structure of the 
formal solution in this case becomes complicated. Two parameters appear 
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and the formal serires solution of the original equation use both of them. 
These two parameters are e as small one and Xk{s) as large one. To simplify 
the representation of the formal series solution in this section we only use 
the segments of the formal asymptotic series. We describe the domain of the 
suitability of this representation of the formal series solution. It is convenient 
to construct the leader term of the intermediate asymptotics in the form 

Ak = Jhi Ak +8'/' Ak ) + ^'/'A,( Ak +e'/' Ak ) + e'^'xT Ak, 



e"'\T{ B, +e'/' B, )• 

The coefficients of the asymptotics contain the independent variable with the 

large term ak = A^' T^, namely, Ak=Ak (c^fe, k, e) and Bk=Bk (c^fc, k, e). 
The leader term of the asymptotics satisfies to the equation 

' 

Ak -2f/* Bk= 0, 
o' „ 

Bk +3t/* Ak = U.- 

It is easy to write the equations for the correction terms of the asymptotics 

Ak -2f/* Bk= 0, 
Bk +6t/* AkAk= U.Tk/Xk] 

lO' 10 

Ak -2f/* Bk= -2U, AkBk, 

lO' 10 5 „ „ „ 

Bk +6C/* AkAk^ -3C/* Ak + Ak - Ak - Bk] 

15' 15 5 5 

Ak -2t/* Bk= -2t/*( AkBk + AkBk ), 

15' 15 , „ 5 Ti, 5 5 10 

Bk +6C/* AkAk^ (1-3 Ak) Ak +T^Ak - 2C/* BkBk - 6C/* AkAk, 

Afc 

20 , 20 „ 

Ak' - 2[/. Bk= - Bk- Ak Bk - 

10 10 1 5 5 , 

2C/*( AkBk + AkBk + ^ AkBk ), 

V '^fe 

20 , 20 10 „ 10 „ 

Bk + 6C/* AkAk=Ak -3 Ak Ak - AkBk - 

10 I , 5 „ 5 „^ 

2C/* BkBk + ^( - 3 AAk - Bk) 
Afe 

1 5 10 2 5 15 

+ -^Tk Ak - 3C/* Ak - 6C/* AkAk ■ 



Afe 
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As was shown above it is convenient to investigate the second order dif- 

n 

ferential equation for the functions Ak- In particular, we have 



,, „ 
Ak +^Ak= 1. 



The solution of this equation matched with the asymptotics of the internal 
formal solution in the neighborhood of the k — 1-th pole is 

jl= -2p(afc,l/2,73(/i;,e)), where -^^{k,e) = 



The parameter g^{k,e) and corrections of the asymptotics are obtained 
as was shown above in the domain of the intermediate expansion. It is easy 
to see that unique type of singularities is connected with the parameter A^. 
The intermediate expansion is valid as 



e"'\k 
vAi, 



< 1, or Afc < £ 



The constructed sequence of the intermediate and the second internal 
expansions is valid in the domain as 



1/6 



6 Fast oscillating asymptotic solution 

In this section we construct a special fast oscillating asymptotic solution of 
the equation @j by Krylov-Bogolyubov method DH, Ijl^. This constructed 
solution is valid as {t^ — t)e~'^/^ ^ 1. This special solution degenerates at the 
point t = t^,. In the neighborhood of this point this special solution loses the 
fast oscillating structure and matchs with the asymptotics from the previous 
section. 



6.1 The family of the fast oscillating solutions 

Here we construct the family of the fast oscillating solutions which depends 
on two parameters. 

The fast oscillating solution of the equation (|^) is seeking in the form 

U{t, e) =U ih, t,e)+eU {ti,t, e) + U {ti,t,e) + . . . , (66) 

where ti = S{t)/e + is the fast variable, the functions S{t) and are 
unknown functions. 
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The problem on constructing the fast oscillating solution for the second 
order differential equation was investigated by many authors. In this section 
we construct the solution in the manner of F.J.Bourland and R.Haberman 
||15|| . They have investigated the second order equation so the direct reference 
for the equation @ is not valid. In this section we represent some evaluations 
which lead us to the elegant formulas from the work [|15] . 



Let us substitute (|66|) in equation (||). As a result we obtain the equations 
for the leader term of the asymptotics ( ]66| ) and the first correction term. 

„ 

iS'dt, U+\U\^U-t U= 1. (67) 

1/0„\10„1 ,0 

zS'dt, U + [2\u\^-t) U + U^U*=-tdtU U . (68) 



2 /,0,„ \2 0„2 

iS'dt, u + h\u \^ -t) u + u*= 

-idt U -i<P'dt, U -2| U?U -{Uf U*. (69) 

The time variable t is the independent variable in these equations. The 
equation (|67|) has the first integral with respect to "fast" variable ti 

1 , „ 

-\U\'-t\u\^-{U + U*)=E{t), (70) 
where E{t) is a "constant" of integration. 

This expression for the first integral can be considered as the equation for 



the function JJ- Let us express U through the complex conjugate function 



and as a result we obtain the equation only for U- 



iS'dt, t/= ±V 2 + {2E{t) +t^)U^ + 2tU +1. 



This equation can be easy integrated 



tS' r , =h + So. (71) 

•^"0 ±^2y3 + (2E + t'^)y^ + 2ty + l 

Where uo and 5*0 are constants and we integrate on the curve r(t) on the 
complex plane y. This curve is defined by 

hy\'-t\y\'-iy + y*) = m- 
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The sign + or — fixed the sheet of a Riemann surface where the initial point 
Uq is situated. We choose the "+" sign. 

The motion defined by (^) is periodic with respect to fast variable ti. 
By integrating of the formula (^) over the period we obtain: 

iS' I = T. (72) 

Mt) /2y3 + (2E + t2)y2 + 2ty + l 



The constant T is the period of the function JJ with respect to fast variable 
ti. This formula is the differential equation for the function S{t). 



It is convenient to represent the leader term JJ (ti,t) as the sum of the 
real and the imaginary parts 



U=Ur +i Ui . 



By shift of phase variable ti we can obtain that Ub. (ti,t) is even function 



with respect to ti and Ui iti,t) is odd one. This property of the functions 



Ui (ti,t) and Ur (ti,t) is convenient to evaluate the integrals with respect 
to fast variable. 



To determine the leader term of the asymptotics ( pq ) it is necessary to 

obtain an equation for E{t). This equation appears from a condition of the 

1 

boundedness of the first correction term JJ of the asymptotics. 

The general solution of the equation for the first correction can be repre- 
sented by two linear independent solutions of the homogeneous equation 

, , ^ OeS „ 

ui[ti,t) = dt^ U, and U2[ti, t) = -^^i^i + Oe U ■ 

In these formulas the Ui is periodic bounded function, but U2 grows linearly 
with respect to ti. The Wronskian of these solutions is 

W = UiU*2 - UIU2 = 

It can be evaluated by explicit formulas 

1 o„o ol 0„0 

W = —{-{\U\^-t)U +l)dE U* + -^{-{\U\^-t)U* + l)dE U= 

1 / 1 , ,4 ,0,9 *N -1 

To obtain necessary condition of the boundedness of the first correction 
terms of the asymptotics (pBD we multiply the equation (RBI) to dt^ JJ * and 
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integrate it with respect to the variable ti over the period. Then we integrate 
the first term by parts 



./■^ 1 f'^/0„ 0^0X1 

iS' dt.dt, U dt, U*^-J^ dh[2{\ U\^- t)dt, U* + ( U*fdt, UjU- 



As a result we obtain: 



Jo 



dti 



0„1 .0,9 01 

U^U*dt, U*-{U*fdt, UU 



fT , fT 

-i / dtidt U dt,U* - / dtidt, U dt,U*. 
Jo Jo 

Let us consider the equation for the complex conjugate function 

1 /,0,„\1 0„1 

-iS'dt, U*+[2\u\^-tju* + { U*f [/= idt U* + i(t>dt, U*. 



Multiply this equation to dt^ JJ and integrate it with respect to ti over the 
period. Then we obtain 



dti 



9 1 „ 1' 

U*) Udu U-U^du U* U* 



fT , l-T 

i / dtidt U*dt, U +1(1)' / dUdt, U dt, U*. 
Jo Jo 

Combining obtained expressions we get 

fT fT 

/ dtidt Udt,U* + i dhdt U*dt, U= 0. 
Jo Jo 

After integrating by parts we get 

fT 



or 



or 



f^ f^ 

; / dti U*dtA U +1 / dhdt U*dt, U= 0, 
Jo Jo 

fT 

idt / dh U*dt, U= 0. 
Jo 

I = i u*du = a = const. 
Jr{t) 



(73) 



The integral in this formula geometrically means the doubled square of the 
domain bounded by the curve r(t). The necessary condition of the bounded- 
ness of the first correction term is invariance of this square. It is the adiabatic 
invariant. 
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The solution of the equation for the first correction term can be written 
in the form 



U (ti, t) = cidt, U +^dE U + R (ti, t). (74) 

Where 

R {ti,t) =ui{ti,t) r dt\dt U {t',t)u;{t',t) + dt U*{t',t)u2{t',t))- 

- U2{tut) r dt'idt U {f, t)ul{t', t) + dt U*{t', t)Ui{t', t)). 
Jo 

To evaluate the second correction we use the property of even and odd of the 

real and imaginary parts of the first correction term. These properties of the 

first and the second terms in the formula ([7^) can be obtained from analysis 

of the leader term of the asymptotics. The last term of this sum contains the 

1 1 
even function ^{R) and the odd function '^{R) with respect to variable ti. 

Constructing of the periodical second correction of the asymptotics leads 
us to determining of the function The necessary condition of this prop- 
erty for the second correction is 

j-T f 2 2 \ 

dt^^Fdt, U*+F*dt,uj =0. 

2 

Where F is the right-hand side of the equation (|69[) . 

Let us separately evaluate the integral depending on square of the term 

1 

U 

fT 1„0 1„0 loO IgO 

J = / dti[ -{2\U?U + U*)dt, U* - (2| U\^U* + { U*f U)dt, U ]. 
Jo 

1 

Integrating by parts and substituting of the expressions for dt^ JJ leads us to 
the following formula 

j = r dt^[{ - idt U -i<P'dt, U)dt,U* + {-idtU* + i<p'dt, U*)dt, U ] 
Jo 

Then we obtain 

fT 10 10 

dt / dti[i U*dt, U -iUdt, U*]. 
Jo 

1 

Substituting of U in this formula and using of even properties gives us 

S' 
dt dhjf-{dEUdt,U*) = 

Jo Oeo 
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-^—^OeI = 01 = const. (75) 

Geo 

This elegant form of the equation for the phase shift was obtained by F. J.Bo- 
urland R.Haberman for the solutions of the second order differential 
equation. 

The formulas (|7TD, ([72D, (^) and ([73| ) determine the family of the fast 
oscillating asymptotic solutions of the equation (j^). The solutions of this 
family contain three parameters. The constant a determines the trajectory 
on the phase plane. The constant 0i in the formula (|75|) and initial condition 
'/'lt=const ~ 'i'o determine the phase shift. The periodic with respect to ti 
leader term of the asymptotics depends on three constants. Specificity of this 
fact is explained by the large values of the time variable. So to determine the 
phase shift correctly we should take into account one condition for the first 
correction term of the constructed asymptotics. The period T is arbitrary 
nonzero constant. 



6.2 The confluent asymptotic solution 

In this section we choose the values of the parameters when the constructing 
fast oscillating solution degenerates at the moment t = t*. The algebraic 
solution constructed above degenerates at the moment t = t^, in the point 
U^. And internal asymptotic solution has as the leader term the constant f/*. 
It leads to the fact that the fast oscillating solution degenerates at the point 
f/* and it allows us to match constructed solutions in the different domains. 
This reason allows us to evaluate the value E{t^:) = E^, 

E* = l\U,\^~t\U,\^ ~{U, + U:). 

By substituting these values we obtain = \ ^ . At the point f/* all 
three roots of the equation 

2|/3 + (2E + t2)y2 ^ 2t,y + 1 = 

coalesce. It means that the integrand in the equation ([7l|) has the singularity 
of the order of 3/2 and elliptic integral in ( [7T| ) degenerates. The obtained 
values of the constants E^ and t^, allow us to evaluate the adiabatic invariant 
a. To obtain this value it is necessary to evaluate the integral 

h = u*du = a^. (76) 
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As result the equation ( |73D defines the function E(t) for our confluent so- 
lution. Thus to obtain the leader term of the asymptotics (^) completely 
it is necessary to determine the phase shift 5*0 in the formula ([7T|). It can 
be evaluated by matching of the asymptotics ([66| ) and internal asymptotics 
from previous section. 



6.3 The domain of suitability of confluent asymptotic 
solution as ^ ^ — 

To match the internal asymptotic expansions with the fast oscillating expan- 
sion it is necessary to know where the last one is valid as t — t*. From the 
formula (^) it is easy to see that the function 5" is defined by the value 
inverse to an integral 



K 



dy 



m ^2y^ + {2E + t^)y^ + 2ty + 1 



So to obtain the domain of suitability of the fast oscillating solution it is 
necessary to determine the order of singularity of the leader term of the 
asymptotics of this integral as t — > — 0. 

The integral K contains two parameters t and E. The function E = E{t) 
is defined by equation (^). In the neighborhood of the point t^^ we use the 
following notation 

t = U+ E = E^ + 6. 
Let us change the variable in the integral 



y 



1\ 1/3 



Then the integral K looks like 

K=[ rfcfe^ + C' 

Jr(t) \ 



26 + 6 



1\ 1/3 



fl + fl 



c 



XX 2/3 



■J 



+ 

+ 

-1/2 



The curve r{t) a.s t = t^: passes through a point C = 0- At this point the 
integrand at the moment t = has nonintegrable singularity and integral 
diverges. 
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The integrand and curve r{t) give us the information that to evaluate 
the main term of the asymptotics as t ^ — of this integral it is enough 
to evaluate the integral only on the small segment of the curve r{t) near the 
point C = 0. 

To make this evaluations it is convenient to pass to the integral on a small 
segment of the real axis and to write the complex variable ( in the form 

C = ^ + ^^, e,^eR. 

In a small neighborhood of point C = the imaginary part of the ( can be 
expressed by real part by formula 

^±(0 = 



± 



21/3 + 22/3^ -e + l^- 2^/^/ 1 + 3 ( ^) + 21/35 + 2-2/3^2 + 24/3^. 



Now if we assume that ^, /i and S are small then we obtain 

•1\ 1/3 



tjI = -2-2/3^3 ^ 2- 



1/3 



^'2 



+ 2^IH 



/i/2 - 2^/35 



+ 



o{e) + + o^e^) + om + o(e/x) + o{i?) + o{5^) + o(/i5). (7?) 

The function stays undetermined in this formula. 

The dependence E{t) and hence can be obtained from (|73|). Let us 
write this equation by integrals on segments 



dx\l t - + + - dx\lt - x2 - Vt^ + 2E + 4x = h. 

ixi Jxi 

in 

Where xi = — (t^ + 2E)/4, x+ is the right real root of the equation 



x^ - 2tx^ - 4a; - 2E = 0, 



(79) 



and X- is the left one. 

To construct the asymptotics of the function E{t) as t t* — it's 
necessary to write the asymptotics of the roots of the equation (^). This 
equation in the investigated domain has two real roots and two complex 
conjugate roots. So we write 



(x — x^)(x — x^)({x — mf' + v?) = 0. 

Where x_, x+, m n are real functions with respect to t. They satisfy to 
Vieta equations 
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X-X+{rn^ + n^) = —2E, 
— (x_ + x+)(m^ + n^) — 2mx_x+ = —4, 
+ + + 2m(x_ + x+) = — 2t, 

x_ + a;+ + 2m = 0. 

We seek a solution of these system of equations in the form 

x_ = ?7* + a, x+ = — 4f/* + /5, m = + rrti. 

The substituting of the values [/*, i?*, t,,, leads us to 

2^/^[3mi - - /i] - 2mi[-mi + + ^] = 0, 
2^/^ [6m? - 2n^ - /i] + 2^/^mi[/i - 3m? - 2n^] + 
+ [m? + n2][-2/i + 3m? - n^] = -26. 

The asymptotics of the solution for this system of the equation is seeking 
in the form 

mi = /iiV^ + /i2/i + . . . , 
n = Uiy/^ + i/2At + . . . , 
S = S^^ + 52{-^f'^ + .... 

By substituting into the system we get 

Thus in the first corrections we did not determine a relations between /xi and 
Vi. This relation can be determine from integral equation (|78D . 

Let us construct the main term of the asymptotics of the integral K. In 
the small neighborhood of the point C = the curve r(t) intersects of the 
real axis. So the integral over a small neighborhood of the point C = can 
be represented in the form 



^" " ^2a-2[-24/3 + 2]K+ + «+(-/i)3/2 

2(3 _ 2[_24/3 + 2]K- + a4-fi)^/^ ' 

where ^rnax is the null of the radicand, C± = ^ + ^^±(0; are constants. 
The parameter a is chosen such that a(— ^ 1 and a -C 1 as /x — > — 0. 
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The asymptotics of the function ri± and the exphcit formula for the main 
term of the integral Iq and change of the variable ^ = (—fiY^'^u in integral 
lead us to the formula 

Ko = 0((-/.)-i/^). 

Therefore 

S' = 0((-/i)^/^), as ^ -0. 
By the similar evaluation one can show that 

OeI = 0((-/x)-i/^), OeS = 0((-/i)i/^), as /i -> -0. 

Therefore 0' = ©((-/i)^/^^. 

By this way we obtain the formula 

dt U= 0(/i~^) asfi^ -0. 

Using Ml and U2 — solutions of the homogeneous equation one can obtain 
that 

As a result we obtain that the fast oscillating asymptotics is valid as 
is easy to see that the domains, where the constructed 
asymptotic solutions are valid, are intersected. 



6.4 Matching of the asymptotics 

Let us match the fast oscillating asymptotics and internal asympotics (O 



(46) into the transition layer. 

The main term of the fast oscillating asymptotics is defined as implicit 
function by the formula 



= -tS' I 

Jv 



dy 



lyem -\y\^y-ty + l 
This integral can be represented in the form as t — + 

t^ = -tS'i / —— ^ — — + 

lr{u) -\y\^y - Uy + 1 + 0{t - Qy 



dy A dy*dy, 

<5Sr 



1 



\y\^y-U + l + 0{t-Qy. 
Where 5Sr is a figure bounded by the curves T(t) and T(t^). 
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The function + W is periodic. In the case the following inequalities 

are 

-C 1 (where 6ar- square of the figure 5Er) < 1, 

we obtain the following representation in some neighborhoods of the point 

t{k) 

Where 

The argument of the function W can be represented in the form 

h_ S{t{k)) + S'{t{k)){t - t(fc)) + o{t - t{k)) (t>{t{k)) . . 

S'^ e{S'{t^k)) + 0{S"{t^k)){t-t^k))) S'i^k))' ^ 

By matching with the main term of the internal asymptotics we obtain 

r 7-- r ~ f)k- 



This matching with the second internal asymptotics gives us that this fast 
oscillating solution also can be matched with the intermediate expansion. 
The constants 0o and 0i evaluate by this way. These constants define the 
phase shift of the fast oscillating solution. But this kind of evaluations are 
not made in this work. 

We are grateful to L.A. Kalyakin, B.I. Suleimanov for stimulating discus- 
sions. And we also thank N. Enikeev for his helpful and comments. 
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